function centralMoment = hwozd_cm( I, p, q )
    momentMi = hwozd_mi(I,p,q);
    centralMoment = momentMi / (I.m00^((p + q + 2)/2));

